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A general technique for the periodic orbit quantization of 
systems with near-integrable to mixed regular-chaotic dynam- 
ics is introduced. A small set of periodic orbits is sufficient 
for the construction of the semiclassical recurrence function 
up to, in principle, infinite length. As in our recent work the 
recurrence signal is inverted by means of a high resolution 
spectral analyzer (harmonic inversion) to obtain the semi- 
classical eigenenergies. The method is demonstrated for the 
hydrogen atom in a magnetic field. To our knowledge this is 
the first successful application of periodic orbit quantization 
in the deep mixed regular-chaotic regime. 

PACS numbers: 05.45. +b, 03.65.Sq, 32.60. +i 



The question of how to semiclassically quantize sys- 
tems with nonintegrable Hamiltonians had puzzled the 
great minds of physics in the early days of quantum me- 
chanics G|. The advent of "exact" quantum mechanics, 
and, later, the availability of more and more powerful 
computing resources, with the possibility of numerically 
solving, in principle, Schrodinger's equation for every 
complex quantum system, pushed that old question in 
the background for several decades. However, the desire 
for deeper physical insight into "what the computer un- 
derstands" Q| finally triggered a renaissance of semiclas- 
sical theory, which persists. It will be the objective of this 
Letter to contribute to solving that longstanding problem 
of semiclassical quantization of nonintegrable systems in 
the mixed regular- chaotic regime. 

The breakthrough of modern semiclassical theory came 
when Gutzwiller proved, by an application of the method 
of stationary phase to the semiclassical approximation 
of the quantum propagator, that for systems with com- 
plete chaotic (hyperbolic) classical dynamics the density 
of states can be expressed as an infinite sum over all (iso- 
lated) periodic orbits, thus laying the foundation of peri- 
odic orbit theory ||. Special methods were designed for 
specific systems to overcome the convergence problems of 
the semiclassical trace formula, i.e., to analytically con- 
tinue its range of convergence to the physical domain 
j|-[7|. None of these methods, however, has succeeded 
so far in correctly describing generic dynamical systems 
with mixed regular-chaotic phase spaces. 

On the other extreme of complete integrability, it is 
well known that the semiclassical energy values can be 
obtained by EBK torus quantization pi. This requires 
the knowledge of all the constants of motion, which are 



not normally given in explicit form, and therefore prac- 
tical EBK quantization based on the direct or indirect 
numerical construction of the constants of motion turns 
out to be a formidable task ||. As an alternative, EBK 
quantization was recast as a sum over all periodic or- 
bits of a given topology on respective tori by Berry and 
Tabor |)j . The Berry- Tabor formula circumvents the nu- 
merical construction of the constants of motion but usu- 
ally suffers from the convergence problems of the infinite 
periodic orbit sum. 

The extension of the Berry- Tabor formula into the 
near-integrable (KAM) regime was outlined by Ozorio de 
Almeida and elaborated, at different levels of refine- 
ment, by Tomsovic et al. and Ullmo et al. pl|| . These 
authors noted that in the near-integrable regime, ac- 
cording to the Poincare-Birkhoff theorem, two periodic 
orbits survive the destruction of a rational torus with 
similar actions, one stable and one hyperbolic unstable, 
and worked out the ensuing modifications of the Berry- 
Tabor formula. In this Letter we go one step further 
by noting that, with increasing perturbation, the stable 
orbit turns into an inverse hyperbolic one representing, 
together with its unstable companion with similar ac- 
tion, a remnant torus. We include the contributions of 
these pairs of inverse hyperbolic and hyperbolic orbits in 
the Berry- Tabor formula and demonstrate for a system 
with mixed regular-chaotic dynamics that this procedure 
yields excellent results even in the deep mixed regular- 
chaotic regime. The system we choose is the hydrogen 
atom in a magnetic field, which is a real physical system 
and has served extensively as a prototype for the inves- 
tigation of "quantum chaos" (for reviews see []L|~[i"4j| ) . 
To our knowledge no semiclassical quantization thus far 
in the mixed regular-chaotic region has previously been 
given. 

The fundamental obstacle bedeviling the semiclassical 
quantization of systems with mixed regular-chaotic dy- 
namics is that the periodic orbits are neither sufficiently 
isolated, as is required for Gutzwiller's trace formula [[|, 
nor are they part of invariant tori, as is necessary for the 
Berry- Tabor formula ||. However, as will become clear 
below, it is the Berry- Tabor formula which lends itself 
in a natural way for an extension of periodic orbit quan- 
tization to mixed systems. We consider systems with a 
scaling property, i.e., where the shape of periodic orbits 
does not depend on the scaling parameter, w = , and 
the classical action S scales as S = sw with s the scaled 
action. For scaling systems with two degrees of freedom, 
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which we will focus on, the Berry- Tabor formula for the 
fluctuating part of the level density reads 
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with M = (Mi,M 2 ) pairs of integers specifying the 
individual periodic orbits on the tori (numbers of ro- 
tations per period, M2/M1 rational), and sm and ?7m 
the scaled action and Maslov index of the periodic or- 
bit M. The function qe in (|l|) is obtained by invert- 
ing the Hamiltonian, expressed in terms of the actions 
(iijib) of the corresponding torus, with respect to I2, 
viz. H{h,h = 9e{I\)) = E The calculation of g% 
from the actions I2) can be rather laborious even for 
integrable and near-integrable systems, and, by defini- 
tion, becomes impossible for mixed systems in the chaotic 
part of the phase space. Here we will adopt the method 
of Ref. 1 11 1 and calculate g' E , for given M = {^1,^2), 
with (fix, [12) coprime integers specifying the primitive 
periodic orbit, directly from the parameters of the two 
periodic orbits (stable (s) and hyperbolic unstable (h)) 
that survive the destruction of the rational torus M, viz. 
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wi4As \ v / det(M s - T) y/- det(M h - I) 



As = -(s h - s s ) 



(2) 
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the difference of the scaled actions, and M s and Mh the 
monodromy matrices of the two orbits. The action sm 
in (llj) is to replaced with the mean action 
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Eq. |^ is an approximation which becomes exact in the 
limit of an integrable system. 

It is a characteristic feature of systems with mixed 
regular-chaotic dynamics that with increasing nonintc- 
grability the stable orbits turn into inverse hyperbolic 
unstable orbits in the chaotic part of the phase space. 
These orbits, although embedded in the fully chaotic part 
of phase space, are remnants of broken tori. It is there- 
fore natural to assume that Eqs. |l| and || can even be 
applied when these pairs of inverse hyperbolic and hy- 
perbolic orbits are taken into account, i.e., more deeply 
in the mixed regular-chaotic regime. 

It should be noted that the difference As between the 
actions of the two orbits is normally still small, and it is 
therefore more appropriate to start from the Berry- Tabor 
formula for semiclassical quantization in that regime than 
from Gutzwiller's trace formula, which assumes well- 
isolated periodic orbits. It is also important to note that 



the Berry- Tabor formula does not require an extensive 
numerical periodic orbit search. The periodic orbit pa- 
rameters S/M2 and g' E are smooth functions of the rota- 
tion number M2/M1, and can be obtained for arbitrary 
periodic orbits with coprime integers (Mi, Ma) by inter- 
polation between "simple" rational numbers M2/M1. 

We now demonstrate the high quality of the exten- 
sion of Eqs. [I] and ^ to pairs of inverse hyperbolic and 
hyperbolic periodic orbits for a physical system that un- 
dergoes a transition from regularity to chaos, namely the 
hydrogen atom in a magnetic field. This is a scaling 



system, with w = 7 
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h cf{ the scaling parameter 



and 7 = B/(2.35 x 10 5 T) the magnetic field strength in 
atomic units. Introducing scaled coordinates 7 2 / 3 r and 
momenta 7 _1 / 3 p and choosing the projection of the an- 
gular momentum on the magnetic field axis L z = one 
arrives at the scaled Hamiltonian 
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with E = Ej~ 2 < 3 the scaled energy. At low energies 
E < —0.6 a Poincare surface of section analysis of the 
classical dynamics [[l3| exhibits two different torus struc- 
tures related to a "rotator" and "vibrator" type motion. 
The separatrix between these tori is destroyed at a scaled 
energy of E w —0.6, and the chaotic region around the 
separatrix grows with increasing energy. At E — —0.127 
the classical phase space becomes completely chaotic. We 
investigate the system at scaled energy E = —0.4, where 
about 40% of the classical phase space volume is chaotic 
(see inset in Fig. [l]), i.e. well in the region of mixed dy- 
namics. We use 8 pairs of periodic orbits to describe 
the rotator type motion in both the regular and chaotic 
region. The results for the periodic orbit parameters 
s/2irM2 and g' E are presented as solid lines in Fig. [l]. The 
squares on the solid lines mark parameters obtained by 
pairs of stable and unstable periodic orbits in the regular 
region of the phase space. The diamonds mark parame- 
ters obtained by pairs of two unstable (inverse hyperbolic 
and hyperbolic) periodic orbits in the chaotic region of 
phase space. The cutoff is related to the winding angle 
4> = 1.278 of the fixed point of the rotator type motion, 
i.e., the orbit perpendicular to the magnetic field axis, 
(M 2 /Mi) cu toff = ir/<f> = 2.458. The solid lines have been 
obtained by spline interpolation of the data points. In 
the same way the periodic orbit parameters for the vi- 
brator type motion have been obtained from 11 pairs of 
periodic orbits (see the dashed lines in Fig. 0). The cutoff 
at M2/M1 = Tr/(j>= 1.158 is related to the winding angle 
<f) = 2.714 of the fixed point of the vibrator type motion, 
i.e., the orbit parallel to the field axis. 

With the data of Fig. | we have all the ingredients 
at hand to calculate the semiclassical density of states 
g(w) in Eq. |l]. The periodic orbit sum includes for both 
the rotator and vibrator type motion the orbits with 
M2/M1 > {M2 /Mi) cutoff- For each orbit the action and 
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FIG. 1. (a) Action s/2irM2 and (b) second derivative 
as a function of the frequency ratio M2/M1 for the rotator 
(solid lines) and vibrator (dashed lines) type motion of the 
hydrogen atom in a magnetic field at scaled energy E — —0.4. 
Inset: Poincare surface of section in semiparabolical coordi- 
nates (/x,p M ; v = 0). 

the function g'^ is obtained from the spline interpola- 
tions. The Maslov indices are tjm = — Mi for the 
rotator and ?7m = 4M2 + 2M% — 1 for the vibrator type 
orbits. However, the problem is to extract the semiclas- 
sical eigenenergies from Eq. [I] because the periodic orbit 
sum does not converge. We therefore adopt the method 
of Refs. |l6|,[l7| where we proposed to adjust the semi- 
classical recurrence signal, i.e., the Fourier transform of 
the weighted density of states w~ 1 ^ 2 g(w) (Eq. [I]) 



C sc (s) = £>W(s-sm) , 
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with the amplitudes being determined exclusively by pe- 
riodic orbit quantities, 
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to the functional form of its quantum mechanical ana- 
logue 



(8) 



where the Wk are the quantum eigenvalues of the scaling 
parameter, and the dk are the multiplicities of the 
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FIG. 2. Semiclassical recurrence signal C sc (s) for the hy- 
drogen atom in a magnetic field at scaled energy E = —0.4. 
Solid and dashed sticks: Signal from the rotator and vibrator 
type motion, respectively. 

eigenvalues (dk — 1 for nondegenerate states). The fre- 
quencies obtained from this procedure are interpreted as 
the semiclassical eigenvalues Wk- The technique used to 
adjust @ to (||) is harmonic inversion p7ft . 

For the hydrogen atom in a magnetic field part of the 
semiclassical recurrence signal C sc (s) at scaled energy 
E = —0.4 is presented in Fig. H. The solid and dashed 
peaks mark the recurrencies of the rotator and vibrator 
type orbits, respectively. Note that C sc (s) can be eas- 
ily calculated even for long periods s with the help of 
the spline interpolation functions in Fig. ^. By contrast, 
the construction of the recurrence signal for Gutzwiller's 
trace formula usually requires an exponentially increas- 
ing effort for the numerical periodic orbit search with 
growing period. 

We have analyzed C sc (s) by the harmonic inversion 
technique in the region < s/2tt < 200. The result- 
ing semiclassical spectrum of the lowest 106 states with 
eigenvalues w < 20 is shown in the upper part of Fig. 
|3|a. For graphical purposes the spectrum is presented as 
a function of the squared scaling parameter w 2 , which is 
equivalent to unfolding the spectrum to constant mean 
level spacing. For comparison the lower part of Fig. ^a 
shows the exact quantum spectrum. The semiclassical 
and quantum spectrum are seen to be in excellent agree- 
ment, and deviations are less than the stick widths for 
nearly all states. The distribution P(d) of the semiclas- 
sical error with d — (w qm — w sc ) / Aw av the error in units 



3 




50 100 150 200 250 300 350 400 

w 2 



T3. 
Q. 




FIG. 3. (a) Semiclassical and quantum mechanical spec- 
trum of the hydrogen atom in a magnetic field at scaled energy 
E — —0.4. (b) Distribution P(d) of the semiclassical error in 



units of the mean level spacing, 
the lowest 106 eigenvalues. 



d = (Wq m - !" S c)/AlB a , for 



of the mean level spacing, Au> av = 1.937/w, is presented 
in Fig. ||b. For most levels the semiclassical error is less 
than 4% of the mean level spacing, which is typical for a 
system with two degrees of freedom [pi] . 

The accuracy of the results presented in Fig. seems 
to be surprising for two reasons. First, we have not ex- 
ploited the mean staircase function N(w), i.e., the num- 
ber of eigenvalues Wk with Wk < w, which is a basic re- 
quirement of some other semiclassical quantization tech- 
niques for bound chaotic systems |^||. Second, as men- 
tioned before, Eq. ^ has been derived for near-integrable 
systems, and is only an approximation, in particular, for 
mixed systems. We have not taken into account any more 
refined extensions of the Berry- Tabor formula (0) as dis- 
cussed, e.g., in Ref. fjj|| . The answer to the second point 
is that the splitting of scaled actions of the periodic orbit 
pairs used in Fig. [j] does not exceed As = 0.022, and 
therefore for states with w < 20 the phase shift between 
the two periodic orbit contributions is wAs = 0.44, at 
most. For small phase shifts the extension of the Berry- 
Tabor formula to near-integrable systems results in a 
damping of the amplitudes of the periodic orbit recur- 
rence signal in Fig. || but seems not to effect the frequen- 
cies, i.e., the semiclassical eigenvalues Wk obtained by the 
harmonic inversion of the function C sc (s). 

In conclusion, we have presented a solution to the fun- 



damental problem of semiclassical quantization of non- 
integrable systems in the mixed regular-chaotic regime. 
We have demonstrated the excellent quality of our pro- 
cedure for the hydrogen atom in a magnetic field at a 
scaled energy E = —0.4, where about 40% of the phase 
space volume is chaotic. The lowest 106 semiclassical and 
quantum eigenenergies have been shown to agree within a 
few percent of the mean level spacings. Obviously, it will 
be straightforward, and rewarding, to apply the method 
to other systems with mixed dynamics. 

This work was supported in part by the Sonder- 
forschungsbereich No. 237 of the Deutsche Forschungs- 
gemeinschaft. J.M. thanks the Deutsche Forschungsge- 
meinschaft for a Habilitandenstipendium (Grant No. Ma 
1639/3). 



[1 
[2 

[3; 

[4; 

[5 

[6 

ir 

[9 
[10 

[11 



[12 
[13. 

[14 



[15 

[is; 
[17; 

[18 



A. Einstein, Verh. Dtsch. Phys. Ges. (Berlin) 19, 82 
(1917); L. Brillouin, J. Phys. Radium 7, 353 (1926); J. 

B. Keller, Ann. Phys. (New York) 4, 180 (1958). 

E. P. Wigner: "It is nice to know that the computer 
understands the problem. But I would like to understand 
it, too." Quoted from E. Heller, Physics Today 46 (7), 
38 (1993). 

M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics (Springer, New York, 1990). 

P. Cvitanovic and B. Eckhardt, Phys. Rev. Lett. 63, 823 
(1989). 

R. Aurich, C. Matthies, M. Sieber, and F. Steiner, Phys. 
Rev. Lett. 68, 1629 (1992). 

M. V. Berry and J. P. Keating, Proc. R. Soc. London 
Ser. A 437, 151 (1992). 

E. B. Bogomolny, Nonlinearity 5, 805 (1992). 
I. C. Percival, Adv. Chem. Phys. 36, 1 (1977). 
M. V. Berry and M. Tabor, Proc. R. Soc. London Ser. A 
349, 101 (1976) and J. Phys. A 10, 371 (1977). 
A. M. Ozorio de Almeida, Hamiltonian Systems: Chaos 
and Quantization (Cambridge University Press, Cam- 
bridge, 1988). 

S. Tomsovic, M. Grinberg, and D. Ullmo, Phys. Rev. 
Lett. 75, 4346 (1995); D. Ullmo, M. Grinberg, and S. 
Tomsovic, Phys. Rev. E 54, 136 (1996). 
H. Friedrich and D. Wintgen, Phys. Rep. 183, 37 (1989). 
H. Hasegawa, M. Robnik, and G. Wunner, Prog. Theor. 
Phys. Suppl. 98, 198 (1989). 

S. Watanabe, Review of Fundamental Processes and Ap- 
plications of Atoms and Ions, ed. C. D. Lin, (World Sci- 
entific, Singapore, 1993). 

O. Bohigas, S. Tomsovic, and D. Ullmo, Phys. Rep. 223, 
43 (1993). 

J. Main, V. A. Mandelshtam, and H. S. Taylor, Phys. 
Rev. Lett. 79, 825 (1997). 

J. Main, V. A. Mandelshtam, G. Wunner, and H. S. Tay- 
lor, Nonlinearity 11, 1015 (1998). 
P. A. Boasman, Nonlinearity 7, 485 (1994). 



4 



